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A Note  on  Improving  Corridor  Flow 
Predictions  in  a Zone  Fire  Model 

Glenn  Forney 


Abstract 

Presently,  the  Consolidated  Fire  and  Smoke  Transport  (CFAST)  zone  fire  model 
assumes  that  smoke  travels  instantly  from  one  side  of  a compartment  to  another.  As 
a result,  upper  layers  start  forming  in  compartments  connected  at  the  end  of  long  cor- 
ridors much  sooner  than  expected.  This  report  documents  a procedure  for  estimating 
ceiling  jet  travel  times  so  that  initial  layer  formation  may  be  predicted  more  accurately 
in  compartments  connected  to  long  corridors.  Standard  correlations  are  not  suitable 
for  predicting  ceiling  jet  velocities  and  temperatures  in  corridors  because  ceiling  jet 
arrival  time  at  a distance  r is  proportional  to  r2  in  ‘normal’  rooms  (rooms  with  length 
to  height  ratios  near  one  ) and  proportional  to  r in  corridors  (rooms  with  large  length 
to  height  ratios).  Correlations  are  derived  by  performing  numerical  experiments  to 
estimate  ceiling  jet  arrival  times,  temperature  fall  off  rates  and  depths  for  cases  with 
various  inlet  ceiling  jet  temperatures  and  depths.  These  correlations  are  used  to  im- 
prove predictions  of  corridor  smoke  flow  in  the  CFAST  zone  fire  model. 


1 Introduction 

A standard  assumption  in  zone  fire  modeling  is  that  once  hot  smoke  enters  a compartment, 
a well  defined  upper  layer  forms  instantly  throughout  the  compartment.  This  assumption 
breaks  down  in  large  compartments  and  long  corridors  due  to  the  time  required  to  fill  these 
spaces.  This  note  proposes  a simple  procedure  for  accounting  for  the  formation  delay  of  an 
upper  layer  in  a long  corridor  by  using  correlations  developed  from  numerical  experiments 
generated  with  the  three  dimensional  Large  Eddy  Simulation  (LES3D)  field  model[l].  Two 
parameters  related  to  corridor  flow  are  then  estimated,  the  time,  tc,  required  for  a ceiling 
jet  to  travel  in  a corridor  and  the  temperature  distribution  down  the  corridor.  These  es- 
timates will  then  be  used  in  the  zone  fire  model,  CFAST,  by  delaying  flow  into  adjacent 
compartments  until  the  ceiling  jet  has  passed  the  compartment. 
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2 Theory 


This  section  outlines  a procedure  for  estimating  a ceiling  jet’s  temperature  decay,  depth, 
velocity  and  hence  arrival  time  at  each  vent  in  a corridor.  A field  model,  LES3D,  is  used  to 
estimate  these  parameters  by  running  a number  of  cases  for  various  inlet  layer  depths  and 
temperatures.  The  vent  flow  algorithm  in  CFAST  then  uses  this  information  to  compute 
mass  and  enthalpy  flow  between  the  corridor  and  adjacent  compartments.  This  is  accom- 
plished by  presenting  the  vent  algorithm  with  a one  layer  environment  (the  lower  layer) 
before  the  ceiling  jet  reaches  the  vent  and  a two  layer  environment  afterwards.  Estimated 
ceiling  jet  temperatures  and  depths  are  used  to  define  upper  layer  properties. 

The  problem  is  to  estimate  the  ceiling  jet  temperature  and  depth  as  a function  of  time 
until  it  reaches  the  end  of  the  corridor.  Two  solution  procedures  were  considered  for  solving 
this  problem.  The  first  involves  the  use  of  the  horizontal  momentum  equation  to  model 
ceiling  jet  flow.  Simpson[2]  gives  a practical  treatment  of  gravity  currents  found  in  nature. 
He  does  not  consider  heat  transfer.  Matsushita  and  Klote[3]  derive  a one  dimensional 
horizontal  momentum  equation,  integrated  over  the  vertical  dimension,  for  tracking  smoke 
flow  in  a corridor.  Rehm  et  al  in  [4,  5]  observes  that  a three  dimensional  model  is  necessary 
for  capturing  the  behavior  induced  by  heat  transfer  between  the  ceiling  jet  and  the  walls. 
Either  approach  is  expensive  computationally  compared  to  the  time  required  to  compute 
other  components  of  a zone  fire  model. 

A second  approach,  the  approach  described  in  this  note,  is  to  run  a field  model  as  a 
pre-processing  step  and  to  summarize  the  results  as  correlations  describing  the  ceiling  jet’s 
temperatures  and  velocities.  An  outline  of  this  process  is  given  by 

1.  Model  corridor  flow  for  a range  of  inlet  ceiling  jet  temperatures  and  depths.  Inlet 
velocities  are  derived  from  the  inlet  temperatures  and  depths. 

2.  For  each  model  run  calculate  average  ceiling  jet  temperature  and  velocity  as  a func- 
tion of  distance  down  the  corridor. 

3.  Correlate  the  temperature  and  velocity  distribution  down  the  hall. 

The  zone  fire  model  then  uses  these  correlations  to  estimate  conditions  in  the  corridor. 
An  outline  of  the  steps  involved  is  given  by 

1.  Estimate  the  inlet  temperature,  depth  and  velocity  of  the  ceiling  jet.  If  the  corridor 
is  the  fire  room  then  use  a standard  correlation.  If  the  source  of  the  ceiling  jet  is 
another  room  then  calculate  the  inlet  ceiling  jet  flow  using  Bernoulli’s  law  for  the 
vent  connecting  the  source  room  and  the  corridor. 

2.  Use  correlations  in  3.  above  to  estimate  the  ceiling  jet  arrival  time  at  each  vent. 
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3.  For  each  vent  in  the  corridor  use  lower  layer  properties  to  compute  vent  flow  before 
the  ceiling  jet  arrives  at  the  vent  and  lower/upper  layer  properties  afterwards. 

2.1  Assumptions 

The  assumptions  made  in  order  to  develop  the  correlations  described  in  this  note  are: 

• The  time  scale  of  interest  is  the  time  required  for  a ceiling  jet  to  traverse  the  length 
of  the  corridor.  For  example,  for  a 100  m corridor  with  1 m/s  flow,  the  characteristic 
time  period  would  be  100  seconds. 

• Cooling  of  the  ceiling  jet  due  to  mixing  with  adjacent  cool  air  is  large  compared  to 
cooling  due  to  heat  loss  to  walls.  Equivalently,  we  assume  that  walls  are  adiabatic. 
This  assumption  is  conservative.  An  adiabatic  corridor  model  predicts  more  severe 
conditions  downstream  in  a corridor  than  a model  that  accounts  for  heat  transfer  to 
walls,  since  cooler  ceiling  jets  travel  slower  and  not  as  far. 

• We  do  not  account  for  the  fact  that  ceiling  jets  that  are  sufficiently  cooled  will  stag- 
nate. Similar  to  the  previous  assumption,  this  assumption  is  conservative  and  results 
in  over  predictions  of  conditions  in  compartments  connected  to  corridors  (since  the 
model  predicts  that  a ceiling  jet  may  arrive  at  a compartment  when  in  fact  it  may 
have  stagnated  before  reaching  it). 

• Ceiling  jet  flow  is  buoyancy  driven  and  behaves  like  a gravity  current.  The  inlet 
velocity  of  the  ceiling  jet  is  related  to  its  temperature  and  depth. 

• Ceiling  jet  flow  lost  to  compartments  adjacent  to  the  corridor  is  not  considered  when 
estimating  ceiling  jet  temperatures  and  depths.  Similarly,  a ceiling  jet  in  a corridor  is 
assumed  to  have  only  one  source. 

• The  temperature  and  velocity  at  the  corridor  inlet  is  constant  in  time. 

• The  corridor  height  and  width  do  not  effect  a ceiling  jet’s  characteristics.  Two  ceil- 
ing jets  with  the  same  inlet  temperature,  depth  and  velocity  behave  the  same  when 
flowing  in  corridors  with  different  widths  or  heights  as  long  as  the  inlet  widths  are 
the  same  fraction  of  the  corridor  width. 

• Flow  entering  the  corridor  enters  at  or  near  the  ceiling.  The  inlet  ceiling  jet  velocity 
is  reduced  from  the  vent  inlet  velocity  by  a factor  of  wvent/wr0om  where  wvent  and 
wroom  are  the  width  of  the  vent  and  room  respectively. 
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2.2  Ceiling  Jet  Characteristics 

2.2.1  Compartments  with  Length  Height  Ratios  Near  One  - Normal  Rooms 

In  a normal  compartment  where  the  length  to  height  ratio  is  near  one,  ceiling  jet  velocities 
can  be  estimated  from  correlations  such  as  ones  found  in  [6]  or  [7].,  or  by  solving  the 
horizontal  momentum  equation  in  addition  to  the  mass  and  energy  conservation  equations. 

Smoke  flow  in  a normal  room  is  qualitatively  different  from  smoke  flow  in  a corridor 
in  one  important  respect.  Corridor  smoke  spreads  in  only  one  dimension,  along  the  length 
of  the  corridor.  Smoke  spreads  in  a normal  room  on  the  other  hand  in  two  dimensions.  In 
addition,  assuming  no  friction  or  heat  transfer  to  boundaries,  ceiling  velocities  in  corridors 
will  be  essentially  constant  while  ceiling  jet  velocities  in  normal  rooms  will  be  approxi- 
mately proportional  to  (1/r) 1 As  a result,  these  correlations  are  not  valid  for  estimating 
ceiling  jet  velocities  in  corridors. 

If  the  length  to  height  ratio  of  the  room  of  fire  origin  is  near  one  then  ceiling  jet  traversal 
times  may  be  estimated  using  the  velocity  correlation  for  a steady  state  fire, 


r/H  < 0.15 
r/H  > 0.15 


(1) 


found  in  [6]  where  H is  the  ceiling  height  (m),  Q is  the  total  energy  release  rate  of  the  fire 
(W)  and  r is  the  distance  from  the  plume  centerline  (m).  The  constants,  0.96  and  0. 195  each 
have  units  of  (m4^3 /(Wl^3s).  CFAST  [8]  estimates  ceiling  jet  temperatures  and  velocities 
using  a correlation  derived  by  Cooper,  given  in  [7]. 

The  time,  t(r),  required  for  the  jet  to  travel  a distance  r from  the  source  can  then  be 
obtained  from  the  velocity  by  integrating  the  quantity  dr/u(r).  Using  the  correlation  given 
in  [6]  we  obtain 


t(r) 


r/H  < 0.15 
r/H  > 0.15 


r/H  < 0.15 
r/H  > 0.15 


(2) 


1 Since  the  surface  area  at  the  exterior  boundary  is  proportional  to  r then  the  velocity  must  be  proportional 
to  1/r  so  that  the  mass  flow  out  remains  equal  to  the  mass  flow  in. 
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Ceiling  jet 

^Ptop  ^Pcj_Pamb^  g<^0 

Tamb 

D • T • 

CJ 

Pamb 

AEw=° 

Figure  1 : Schematic  of  a gravity  current  defining  terms  used  to  estimate  its  inlet  velocity. 

where  f = r/(.lbH)  and  Q = [Q / H)l^z / H . The  time  t(r)  corresponding  to  Cooper’s 
velocity  correlation  must  be  integrated  numerically.2  The  arrival  time  is  approximately 
proportional  to  r2,  the  distance  squared.  It  will  be  shown  later  empirically  that  in  a corridor 
the  arrival  time  is  proportional  to  r,  the  distance. 

2.2.2  Compartments  with  Large  Length  Height  Ratios  - Corridors 

The  correlations  defined  by  (1)  and  (2)  are  not  appropriate  in  corridors.  Ceiling  jet  flow  in  a 
corridor  can  be  characterized  as  a one  dimensional  gravity  current.  To  a first  approximation, 
the  velocity  of  the  current  depends  on  the  difference  of  the  density  of  the  gas  located  at  the 
leading  edge  of  the  current  and  the  gas  in  the  adjacent  ambient  air.  The  velocity  also 
depends  on  the  depth  of  the  current  below  the  ceiling.  A simple  formula  for  the  gravity 
current  velocity  may  be  derived  by  equating  the  potential  energy  of  the  current,  mgd0/ 2, 
measured  at  the  half-height  do/ 2 with  its  kinetic  energy,  mU2/ 2 to  obtain 

U = ^do  (3) 

where  m is  mass,  g is  the  acceleration  of  gravity,  H is  the  height  of  the  gravity  current  and 
U is  the  velocity.  When  the  density  difference,  A p = pamb  ~ Pcj , between  the  current  and 

the  ambient  fluid  is  small,  the  velocity  U is  proportional  to  ^Jgdoj2  = yj gd0  where 

pamb , Tamh  are  the  ambient  density  and  temperature,  pCj,  TCJ  are  the  density  and  temperature 
of  the  ceiling  jet,  AT  = TCj  — Tamb  is  the  temperature  difference.  Here  use  has  been  made 
of  the  ideal  gas  law,  Pamb^amb  ~ PcjTcj. 

This  can  be  shown  using  terms  defined  in  Figure  1 by  using  an  integrated  form  of 
Bernoulli’s  law  noting  that  the  pressure  drop  at  the  bottom  of  the  ceiling  jet  is  Pb  — 0,  the 
pressure  drop  at  the  top  is  Pt  = gdo{pCJ  — Pamb ) and  using  a vent  coefficient  cvent  = 0.74, 

2There  is  no  analytic  formula  for  t(r)  since  the  calculation  depends  on  a correlation  that  is  obtained 
numerically. 
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to  obtain 


U0 


Cvent 


Cvent 


Cyent 
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V8 
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\/8 
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1 Pt  + yj  PfPb  + Pb 

y/Pcj  VPt  + \[Pb 
yPt/pcj 


Pamb  Pcj 
Pcj 


0.7 


AT 

Pamb 


(4) 


Formulas  of  the  form  of  equation  (4)  lead  one  to  conclude  that  a ceiling  jet’s  char- 
acteristics in  a corridor  depend  on  on  its  depth,  do,  and  relative  temperature  difference, 
AT / Tamb • Therefore,  as  the  jet  cools,  it  slows  down,  due  to  the  factor  AT/Tarn&.  If  no  heat 
transfer  occurs  between  the  ceiling  jet  and  the  surrounding  walls,  then  the  only  mechanism 
for  cooling  is  mixing  with  surrounding  cool  air. 


2.3  Numerical  Experiments 

2.3.1  The  Parameter  Study 

In  order  to  better  understand  the  effects  of  the  inlet  ceiling  jet  temperature  and  depth  on 
its  characteristics  downstream  in  a corridor,  a number  of  numerical  experiments  were  per- 
formed using  the  field  model  LES3D[1].  Twenty  cases  were  run  with  five  different  inlet 
depths  and  four  different  inlet  temperatures. 

Boundary  Conditions  - corridor  flow  entrance  The  inlet  ceiling  jet  temperature  rise, 
AT0,  and  depth,  do,  were  used  to  define  an  inlet  velocity,  Uq  using  Uo  = 0.7ygdoAT0/Tam&. 
This  formula  is  derived  from  Bernoulli’s  law  for  flow  through  a vent. 

The  inlet  ceiling  jet  depths,  do,  used  in  the  parameter  study  are  0.15  m,  0.30  m,  0.45 
m,  0.60  m and  0.75  m.  The  inlet  ceiling  jet  temperature  rises,  AT0,  used  in  the  parameter 
study  are  100  °C,  200  °C,  300  °C  and  400  °C.  Velocities  using  equation  (4)  corresponding 
to  these  inlet  depths  and  temperature  rises  are  given  in  Table  1 . 

Boundary  Conditions  - all  other  surfaces  For  each  case,  a non-slip  velocity  boundary 
condition  was  imposed  at  all  solid  boundaries.  Adiabatic  thermal  boundary  conditions  were 
imposed  at  the  walls  to  simulate  no  heat  transfer  to  wall  surfaces.  A vertical  symmetry 
plane  along  the  centerline  of  the  corridor  was  used  to  reduce  the  number  of  grids,  thereby 
improving  the  resolution.  An  open  boundary  condition  was  imposed  at  the  far  end  of  the 
corridor. 
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Table  1 : Inlet  Velocity  as  a function  of  Inlet  Temperature  and  Depth 


Depth  (m) 

Temperature  Excess  (°C) 
+100  +200  +300  +400 

0.15 

0.50 

0.71 

0.87 

1.00 

0.30 

0.71 

1.00 

1.22 

1.41 

0.45 

0.87 

1.22 

1.50 

1.73 

0.60 

1.00 

1.41 

1.73 

2.00 

0.75 

1.12 

1.58 

1.94 

2.24 

Geometry  The  simulated  corridor  had  dimensions  of  10  m x 2.4  m x 2.4  m.  Each  grid 
cell  had  dimensions  of  approximately  10  cm  x 5 cm  x 2.5  cm  (4.0  in  x 2.0  in  x 1.0  in) 
where  the  longest  grid  dimension  occurred  along  the  length  of  the  corridor.  Approximately 
220,000  grid  cells  were  used  to  model  the  corridor;  24  across  the  width  of  the  corridor  and 
96  along  both  the  height  and  length  of  the  corridor.  For  most  cases,  a Reynold’s  number  of 
9200  (approximately  962)  was  used  to  resolve  the  small  scale  flow  features. 

A coarser  grid  with  dimensions  10  cm  x 5 cm  x 5 cm  was  used  initially.  It  was  found 
that  the  vertical  grid  dimension,  5 cm  was  not  sufficiently  small  to  resolve  the  thin,  0.15m 
ceiling  jet  cases. 


Qualitative  numerical  corridor  flow  experiments  Figures  2 through  7 are  presented  to 
illustrate  how  various  boundary  conditions  (adiabatic  walls,  cold  walls  etc ) affect  qualita- 
tive ceiling  jet  characteristics.  These  simulations  were  performed  using  LES3D[1].  Next  a 
parameter  study  was  performed  in  order  to  study  the  quantitative  behavior  of  the  ceiling  jet 
by  varying  the  inlet  ceiling  jet  depth  and  temperature  rise. 

As  can  be  seen  from  Figures  2 and  3 a ceiling  jet  modeled  in  an  enclosure  with  adiabatic 
walls  travels  farther  and  faster  than  a ceiling  jet  modeled  with  cold  walls. 

Figure  4 shows  shaded  temperature  contours  in  a vertical  plane  along  the  centerline  of 
a 10  m corridor  at  10  and  20  seconds.  The  temperature  distribution  at  these  two  times  is 
about  the  same.  Even  though  an  individual  portion  of  smoke  has  moved  up  to  10  m,  the 
peak  temperature  (edge  of  darkest  contour)  has  only  shifted  by  about  one  meter,  from  5 m 
to  4 m. 

Similarly,  Figure  5 shows  shaded  velocity  contours  in  a vertical  plane  along  the  cen- 
terline of  a 10  m corridor  at  20  and  40  seconds.  The  velocity  distribution  does  not  change 
much  between  the  two  times.  Because  of  this,  we  assume  that  the  walls  do  not  heat  up  and 
that  the  temperature  and  velocity  distributions  quickly  reach  steady  state. 

Figure  6 shows  shaded  velocity  contours  in  a vertical  plane  along  the  centerline  of  a 10 
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(b)  cold  walls 
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Figure  2:  Shaded  temperature  contours  along  the  corridor  centerline  for  three  temperature 
boundary  conditions.  The  inlet  flow  has  a depth  of  0.6  m (2  ft),  a velocity  of  1.0  m/s  and  a 
temperature  rise  of  300  °C  above  ambient.  A no-slip  boundary  condition  is  imposed  at  all 
wall  surfaces. 
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(a)  free-slip,  adiabatic 
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(b)  no-slip,  adiabatic 


SUCE  23  IN  XZ  PLANE  FOR  CASE  HALL2N.C0L0. NOSUP 


(c)  no-slip,  cold  walls 


Figure  3:  Shaded  U velocity  contours  along  the  corridor  centerline  for  three  different 
boundary  conditions.  The  inlet  flow  has  a depth  of  0.6  m (2  ft),  a velocity  of  1.0  m/s 
and  a temperature  rise  of  300  °C  above  ambient. 
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Figure  4:  Shaded  temperature  contours  along  the  corridor  centerline  at  two  different  times 
during  the  simulation.  The  inlet  flow  has  a depth  of  0.6  m (2  ft),  a velocity  of  1.0  m/s  and 
a temperature  rise  of  300  °C  above  ambient. 
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Figure  5:  Shaded  U velocity  contours  along  the  corridor  centerline  at  two  different  times 
during  the  simulation.  The  inlet  flow  has  a depth  of  0.6  m (2  ft),  a velocity  of  1.0  m/s  and 
a temperature  rise  of  312  °C  above  ambient. 
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m corridor  for  three  inlet  ceiling  jet  depths.  The  inlet  velocity  in  each  case  is  1.0  m/s.  Note 
that  the  ceiling  jet  velocity  increases  as  the  inlet  depth  increases. 

Figure  7 shows  shaded  velocity  contours  in  a vertical  plane  along  the  centerline  of  a 
10  m corridor  for  three  inlet  ceiling  jet  temperatures.  Note  that  the  ceiling  jet  velocity 
increases  as  the  inlet  temperature  increases. 

2.3.2  Reducing  the  Data 

Estimating  Zone  Fire  Modeling  Quantities  from  Field  Modeling  Results  The  com- 
putational fluid  dynamic  or  field  model,  LES3D,  calculates  temperatures,  pressures  and 
velocities  at  many  points  in  a 3-D  rectangular  grid.  The  zone  fire  model,  CFAST,  on  the 
other  hand  calculates  pressure  at  the  floor,  temperatures  in  the  lower  and  upper  layer  and  a 
layer  interface  height.  In  this  section,  a procedure  is  presented  for  reducing  the  data  from 
a field  model  to  a form  consistent  with  a zone  fire  model.  For  each  vertical  plane  along  the 
length  of  a corridor,  a layer  interface  height  is  estimated  by  using  the  distance  above  the 
floor  where  the  temperature  gradient  is  greatest.  An  upper  layer  temperature  is  calculated 
by  averaging  all  temperatures  in  the  slice  above  the  estimated  layer  height. 

The  simulation  of  these  cases  resulted  in  an  extensive  set  of  velocity,  temperature  and 
pressure  predictions.  This  data  is  reduced  to  a more  manageable  size  by  noting  that  the 
inlet  high-temperature  ceiling  jet  flow  stratifies  the  corridor  gases  into  two  regions;  an 
upper  region  of  hot,  fast  flowing  air  and  a lower  region  of  relatively  cool  quiescent  air.  An 
average  upper  layer  temperature  and  layer  height  is  estimated  for  each  vertical  slice  along 
the  length  of  the  corridor.  To  estimate  the  upper  layer  temperature,  Tavg(i),  and  interface 
height,  yiay(i)  for  each  vertical  slice,  z,  downstream  in  a corridor,  perform  the  following 
calculations. 

The  terms,  ibar,  jbar  and  kbar  refer  to  the  number  of  grid  cells  in  the  i,  j and  k direc- 
tions respectively.  The  z direction  is  down  the  longest  side  of  the  corridor,  the  j direction  is 
across  the  corridor  and  the  k direction  is  along  the  vertical  direction.  Let  T(i,j,  k)  be  the 
temperature  of  the  grid  cell  (z,  j,  k).  Consider  the  z’th  vertical  slice.  For  each  vertical  ele- 
vation k in  the  z’th  slice,  define  an  average  temperature  Tavg , a temperature  gradient,  ATavg 
(ratio  of  the  difference  between  adjacent  average  temperatures  in  the  vertical  direction  and 
the  grid  spacing).  Finally  define  a maximum  temperature  gradient  A Tmax.  Calculate  these 
terms  using 


Tavg  (L 

ATavg{i , k) 
ATmax(i) 


1 

jbar 


jbar 


3= 1 


Tavg{i,k  + 1)  - Tavg{i,  k) 

max  A7^t2Q2;(z,  /c) 
i</t<kbar 


(5) 
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(c)  depth  = 0.60  m 


Figure  6:  Shaded  U velocity  contours  along  the  corridor  centerline  for  three  different  inlet 
ceiling  jet  depths.  The  inlet  flow  has  a velocity  of  1.0  m/s  a temperature  rise  of  312  °C 
above  ambient  and  a depth  of  0.30  m,  0.45  m or  0.60  m 
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Figure  7 : Shaded  U velocity  contours  along  the  corridor  centerline  for  three  different  inlet 
ceiling  jet  temperature  rises,  ATo-  The  inlet  flow  has  a depth  of  0.60  m and  a velocity  of 
1.0  m/s  . 
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ATavg(i,3)=Tavg(i,4)-Tavg(i,3) 
ATavg(i,2)=Tavg(i,3)-Tavg(i,2) 
ATavg(i.  1 )=Tavg(i,2)-Tavg(i,  1 ) 


ATmax(i)  = max(ATavg(i,l),ATavg(i,2),ATavg(i,3)) 

Figure  8:  Illustration  of  an  average  temperature  calculation  for  one  4x4  vertical  slice.  The 
estimated  layer  interface  is  the  row  where  the  greatest  temperature  change  occurs. 


The  layer  height  is  then  the  index  kmax  that  corresponds  to  the  maximum  temperature 
change  given  by  A Tmax.  In  other  words,  the  height  above  the  floor  where  the  temperature 
changes  from  a cool  lower  layer  to  a warm  upper  layer.  The  average  upper  layer  tempera- 
ture for  slice  i is  then  given  by 


Tu(i)  = 


jbar  kbar 

E E T(u,k) 


jbar  (kbar  + 1 - kmal)  jrj  k 
kbar 

E Tavg(i,  k)  . 


kbar  + 1 — k 


max  k=kT 


Figure  8 gives  an  example  of  this  calculation  for  JBAR  = KBAR  = 4. 

The  temperature  rise  above  ambient  for  each  slice  i is  given  by  A T(i)  = Tu(i)  — Tamb. 
These  temperature  rises  are  scaled  by  the  inlet  temperature  rise,  A T0,  and  transformed 
using  log(Ar/AT0).  The  resulting  data  are  presented  in  Figure  9.  Note  that  each  plot  is 
nearly  linear  and  that  all  plots  (except  for  the  d0  = 0.15m  group)  lie  within  a single  group. 
This  implies  that  the  relative  temperature  falloff  is  independent  of  the  inlet  temperature  rise 
and  depth  (assuming  that  the  inlet  depth  is  sufficiently  thick). 

The  temperature  curves  presented  in  Figure  9 were  approximated  by  straight  lines  using 
a linear  least  squares  curve  fitting  procedure. 
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2.3.3  Summary  of  Results 


Ceiling  Jet  Temperatures  The  line  that  fits  the  temperature  fall  off  data  in  a least  squares 
sense  is  given  for  the  plots  in  Figure  9.  This  fit  is  given  in  the  form  of 


This  is  equivalent  to 


where  C\  = 10a  and  hi/2  = — log(2)/6.  The  parameter  h\j2  in  equation  (8)  has  a physical 
interpretation.  It  is  the  distance  down  the  corridor  where  the  temperature  rise  AT,  falls  off 
to  50  per  cent  of  its  original  value  or  equivalently,  A T(x  + /lo.s)  = 0.5AT(x).  Equation 
(7)  and  (8)  are  equivalent  mathematically.  Equation  (8)  is  more  intuitive  however,  because 
the  coefficient,  hi/2  has  a physical  interpretation. 

The  half-distance,  hi/2,  can  be  approximated  by  h\j2  = log(2)/0.018  ^ 37  where 
b = —0.018  is  given  in  Figure  9.  Similarly,  the  coefficients  C\  is  approximated  by  C\  = 
10G  = 10-°-003  ^ 1 where  a is  also  given  in  Figure  9.  Therefore  the  temperature  rise,  AT, 
may  be  approximated  by 


AT  = AT0 


(9) 


Ceiling  Jet  Arrival  Times  Numerical  thermocouples  were  placed  0.15  m (6  in)  below 
the  ceiling  every  0.5  m along  the  center  line  of  the  corridor  for  all  cases  except  for  the 
d0  = 0.15m  cases  where  they  were  placed  0.075  m below  the  ceiling.  Ceiling  jet  arrival 
times  were  recorded  by  noting  when  the  temperatures  rose  1°C.  The  arrival  time  for  each 
case  was  scaled  by  the  ceiling  jet  velocity,  Uo,  at  the  left  entry  vent.  These  reduced  data  are 
displayed  in  Figure  10.  Relative  temperature  fall  off  data  for  all  ceiling  jets  are  displayed 
in  Figure  9.  Note  that  most  arrival  time  curves  lie  within  approximately  the  same  region. 
A group  of  curves,  corresponding  to  the  inlet  depth  of  d0  = 0.15m,  are  separate  from  the 
main  group.  This  is  because  the  0.15m  ceiling  jets  are  weaker  than  the  rest.  They  lose 
their  driving  potential  resulting  in  lower  velocities  and  hence  greater  arrival  times.  The 
arrival  time  of  the  ceiling  jet  head  may  be  approximated  by  t — x/Uo  for  all  d0  except  for 
do  = 0.15m.  As  stated  before,  the  do  = 0.15m  ceiling  jets  are  not  strong  enough  to  sustain 
a well  defined  flow  for  the  length  of  the  corridor. 
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Figure  9:  Common  log  of  relative  temperature  excess  downstream  in  a corridor  using  an 
adiabatic  temperature  boundary  conditions  for  several  inlet  depths  and  inlet  temperature 
rises.  The  inlet  velocity,  U0,  is  given  by  equation  (4) 


3 Summary  and  Conclusions 


A procedure  for  estimating  ceiling  jet  travel  times  using  correlations  has  been  described. 
The  procedure  allows  initial  layer  formation  to  be  predicted  more  accurately  in  compart- 
ments connected  to  long  corridors.  These  correlations  were  derived  by  performing  numer- 
ical experiments  using  the  computational  field  model  LES3D  to  estimate  ceiling  jet  arrival 
times,  temperature  fall  off  rates  and  ceiling  jet  depths  for  cases  with  various  specified  inlet 
ceiling  jet  temperatures  and  depths.  These  correlations  were  then  used  to  improve  predic- 
tions of  corridor  smoke  flow  in  the  zone  fire  model,  CFAST 

The  numerical  experiments  with  LES3D  demonstrated  that  for  the  cases  simulated,  ceil- 
ing jet  characteristics  depend  on  the  relative  inlet  temperature  rise  and  not  the  inlet  depth. 
Flow  in  long  corridors  ( greater  than  10  m)  need  to  be  better  characterized  due  to  the  flow 
stagnation  which  may  occur  because  of  the  ceiling  jet’s  temperature  decay.  Even  though  the 
use  of  LES3D  served  as  a good  tool  to  design  a corridor  flow  sub-model,  comparisons  still 
need  to  be  made  between  real-scale  experiments  and  the  zone  fire  model  predictions.  Such 
comparisons  will  allow  us  to  assess  the  validity  of  the  modeling  assumptions,  to  determine 
the  accuracy  of  the  predictions  and  to  gain  better  confidence  in  the  use  of  this  model. 
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Figure  10:  Ceiling  jet  arrival  time  as  a function  of  distance  down  the  corridor  scaled  by  the 
initial  ceiling  jet  velocity,  Uo  for  several  initial  inlet  depths,  do,  and  temperature  rises,  ATo. 
The  walls  are  adiabatic.  The  inlet  velocity,  Uo,  is  given  by  equation  (4).  The  ceiling  jet’s 
arrival  time  is  measured  by  noting  when  its  temperature  rises  1 °C  above  ambient  at  a given 
distance  downstream  from  the  inlet. 
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A Implementation 

A.l  Data  Structures 

Several  new  variables  have  been  introduced  to  record  the  status  of  the  corridor  model. 
IZHALL  indicates  whether  a room  is  a hall  and  whether  user  defined  velocities  and  depths 
should  be  used.  IZVENT  records  distance  information  for  vents  connected  to  halls.  ZZHALL 
records  times,  velocities  and  depths  ceiling  jets  in  each  corridor.  More  details  are  given  be- 
low. 


HALLDIST(I,J,K,L)  First  note  that  the  triple  (I,J,K)  defines  a unique  vent  in  CFAST. 

HALLDIST(I,J,K,1)  is  then  the  distance  between  vent  (I,J,K)  and  a ref- 
erence point  in  room  I.  It  doesn’t  matter  where  this  point  is  as  long  as 
each  vent  connected  to  ROOM  I uses  the  same  reference  point.  HALLD- 
IST(I,J,K,2)  is  the  distance  between  vent  (I,J,K)  and  a reference  point  in 
room  J.  Similarly,  the  reference  point  for  ROOM  J must  be  the  same  for 
all  vents  connected  to  ROOM  J.  This  variable  is  set  in  NPUTQ. 

IZHALL(IROOM,index)  This  array  records  integer  information  about  a hall  in  room 
IROOM.  The  parameter,  index  may  be  one  of  IHDEPTHFLAG,  IH- 
HALLFLAG,  IHMODE,  IHROOM,  IHVELFLAG,  IHVENTNUM  or 
IHXY  which  are  defined  below. 


IHDEPTHFLAG  If  IZH  ALL(I,IHDEPTHFL  AG)=  1 then  use  the  depth  specified  on 

the  HALL  command  line  for  the  ceiling  jet  depth.  This  variable  is 
set  in  NPUTQ. 

IHHALFFLAG  If  IZH  ALL(I,IHHALFFLAG)  is  non-zero  then  use  the  temperature 

decay  parameter  specified  on  the  HALL  command  line.  This  vari- 
able is  set  in  NPUTQ. 


IHMODE  The  variable  IZHALL(I, IHMODE)  can  have  a value  of  IHBEFORE, 

IHDURING  or  IH AFTER.  If  IHBEFORE,  then  flow  has  not  started 
in  corridor  I.  If  IHDURING,  then  flow  has  started  but  has  not  reached 
the  end  of  the  corridor.  If  IHAETER,  then  flow  has  reached  the  end 
of  the  corridor.  This  variable  is  set  in  SETHALL,  ENITMM  and 
NPUTQ. 

IHROOM  If  IZH ALL(I,IHROOM)=  1 then  room  I is  a hall.  This  variable  is  set 

in  NPUTQ. 
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IHVELFLAG  If  IZH  ALL(I,IHVELFL  AG)=  1 then  use  the  velocity  specified  on  the 

HALL  command  line  for  the  ceiling  jet  velocity.  This  variable  is  set 
in  NPUTQ. 

IHVENTNUM  If  IZHALL(I JHVENTNUM)  is  non-zero  then  it  contains  the  num- 


IHXY 

ber  of  the  vent  that  initiated  flow  into  the  hall.  This  variable  is  set  in 
NPUTQ. 

If  IZHALL(I,IHXY)  is  one,  then  the  X direction  corresponds  to  the 
longest  dimension  of  hall  I.  If  this  value  is  two  then  the  Y direction 
corresponds  to  the  longest  dimension  of  hall  I.  This  variable  is  set  in 
ENITAMB. 

IZVENT(IV,4) 

If  IZVENT(IV,4)=1  then  the  ‘from’  room  is  a hall.  The  from  room  is 
the  first  room  specified  on  the  HVENT  command  line.  This  variable  is 
set  in  DATACOPY. 

IZ  VENT (TV, 5) 

If  IZVENT(IV,5)=1  then  the  ‘to’  room  is  a hall.  The  ‘to’  room  is  the 
second  room  specified  on  the  HVENT  command  line.  This  variable  is 
set  in  DATACOPY. 

UPDATEHALL 

If  UPDATEHALL  is  TRUE  then  the  hall  flow  velocity,  temperature 
and  depth  is  updated  in  the  subroutine  HFLOW.  The  variable  UPDATE- 

HALL  is  set  in  the  subroutine  SOLVE  after  RESID  is  called. 
ZZHALL(I,index)  This  array  records  floating  point  information  about  a hall  in  room 


IROOM.  The  parameter,  index,  may  be  one  of  IHDEPTH,  IHDIST,  IH- 
HALF, IHMAXLEN,  IHORG,  IHTEMP  IHTIMEO  or  IHVEL,  which 
are  defined  below. 

IHDEPTH 

If  IZHALL(IJHDEPTHFLAG)  is  not  zero,  then  the  variable 
ZZHALL(I, IHDEPTH)  contains  the  corridor  flow  depth  specified  on 
the  HALL  command  line  for  room  I.  If  IZHALL(I,IHDEPTHFLAG) 
is  zero,  then  ZZHALL(IHDEPTH)  contains  the  corridor  flow  depth 
calculated  by  the  hybrid  model. 

IHDIST 

The  variable,  ZZHALL(IJHDIST),  contains  the  distance  that  the 
ceiling  jet  has  traveled  down  the  hall.  This  variable  is  set  in  SETHALL. 

IHHALF 

If  IZHALL(I,IHHALFFLAG)  is  not  zero  then  ZZHALL(LIHHALF) 
contains  the  distance  where  the  ceiling  jet  flow  temperature  is  re- 
duced 50  per  cent.  If  IZHALL(I,IHHALFFLAG)  is  zero  then  the 
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variable  ZZHALL(I,IHHALF)  contains  the  50  per  cent  distance  pa- 
rameter specified  on  the  HALL  command  for  room  I.  This  variable 
is  set  in  NPUTQ  and  HALLTRV. 

IHMAXLEN  This  variable  contains  the  maximum  length  that  flow  can  travel  in 


IHORG 

room  I before  it  encounters  the  end  of  the  corridor.  This  variable  is 
set  in  SETHALL. 

The  variable,  ZZ(I, IHORG),  contains  the  distance  down  the  hall 
where  flow  starts.  This  dimension  corresponds  to  either  a vent  open- 
ing or  a fire  location,  where  corridor  flow  originates.  Note  that 
the  variables  ZZHALL(I,IHORG)  and  ZZHALL(I, IHMAXLEN)  al- 
ways sum  to  the  length  of  the  corridor.  This  variable  is  set  in  SETHALL. 

IHTEMP 

The  variable  ZZHALL(I, IHTEMP)  contains  the  initial  temperature 
of  the  ceiling  jet  flow  in  room  I.  This  variable  is  set  in  SETHALL. 

IHTIMEO 

If  ZZHALL(I,IHTIME0)>0  then  this  variable  contains  the  time  when 
smoke  started  flowing  into  hall  I.  This  variable  is  set  in  SETHALL. 

IHVEL 

If  IZHALL(I,IHVELFLAG)  is  not  zero,  then  ZZHALL(IHVEL)  con- 
tains the  corridor  flow  velocity  specified  on  the  HALL  command  line 
for  room  I.  If  IZHALL(I,IHVELFLAG)  is,  then  ZZHALL(IHVEL) 
contains  the  corridor  flow  velocity  calculated  by  the  hybrid  model. 
This  variable  is  set  in  NPUTQ  and  SETHALL. 

ZZVENTDIST(IJ)  The  variable  ZZVENTDIST(I,J)  (if  greater  than  zero)  contains  the 
distance  between  vent  J and  the  vent  where  smoke  started  entering  HALL 
I.  Note  that  vent  J must  be  connected  to  room  I.  This  variable  is  set  in 
SETHALL. 

A.2  Routines 

Three  new  routines  have  been  created  to  update  hall  data  structures  and  to  calculate  corre- 
lations for  ceiling  properties  in  halls. 


HALLHT 

This  routine  computes  the  temperature  and  velocity  of  the  ceiling  jet  at 
each  detector  location  in  a corridor.  The  routine  CEILHT  performs  the 
same  function  for  rooms  that  are  not  corridors.  The  correlations  are  used 
only  while  the  ceiling  jet  is  traveling  down  the  corridor. 
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Figure  1 1 : Diagram  for  the  HALLDIST  command 


HALLTRV  This  routine  computes  the  ceiling  temperature,  density  and  velocity  at 
a given  distance  from  where  smoke  enters  the  corridor.  The  routine 
HALLHT  calls  HALLTRV  for  each  detector  location. 

SETHALL  This  routine  updates  the  hall  data  structures,  ZZHALL  and  IZHALL 
each  time  control  is  returned  from  DASSL.  In  particular,  it  re-calculates 
ceiling  jet  velocities,  temperatures  and  determines  whether  the  ceiling 
jet  has  reached  the  end  of  the  corridor. 


A.3  Algorithms 

This  section  presents  an  outline  of  the  calculations  that  were  performed  to  implement  the 
hybrid  model  by  relating  the  equations  and  correlations  developed  in  previous  sections  with 
FORTRAN  code  used  to  perform  the  calculations. 

1 . Initialize  various  components  of  ZZHALL  and  IZHALL  in  the  subroutine  INITMM. 

2.  Read  in  the  HALL  and  HVENT  commands  to  define  which  rooms  are  corridors  and 
where  vents  are  located  in  the  hall. 


22 


3.  Update  the  hall  data  structures.  Hall  data  structures  are  updated  after  DASSL  com- 
pletes a time  step.  In  SOLVE,  the  logical  variable  UPDATEHALL  is  set  allowing 
SETHALL  to  be  called  from  both  the  subroutines,  FIRES  and  HFLOW  (a  ceiling  jet 
in  a hall  may  start  due  to  a fire  in  the  hall  or  flow  through  a vent  from  an  adjacent 
room). 

(a)  In  FIRES,  determine  whether  the  room  containing  the  main  fire  is  a hall.  If  it 
is  a hall,  then  use  the  ceiling  jet  routine  CEILHT,  to  estimate  initial  ceiling  jet 
temperature,  velocity  and  depth. 

(b)  In  HFLOW,  determine  whether  smoke  has  started  flowing  from  an  adjacent 
room  to  a hall.  If  it  has,  then  record  the  start  time  and  vent  number  into  the 
ZZHALL  and  IZHALL  data  structures.  HALL  is  flowing  into  the 

A.4  New  Input  Command 

In  order  to  specify  CFAST  cases  involving  corridor  flow  a new  input  command  HALL  was 
created  and  the  command  HVENT  was  modified.  The  command  used  to  define  a room  as 
a hall  is 

HALL  iroom[,  vel,  depth,  tdecay] 

where  iroom  is  the  room  number,  and  vel,  depth  and  tdecay  are  optional  parameters 
used  to  specify  the  ceiling  jet  velocity,  depth  and  distance  where  the  temperature  falls  off 
by  50  per  cent.  The  HVENT  input  command  was  modified  to  allow  the  user  to  enter  further 
dimensional  data  about  a vent. 

The  new  HVENT  command  is  given  by 
HVENT  I,  J,  K,  W,  H,  L,  VWCOS,  HI,  HJ 

The  parameters  I,  J,  K,  W,  H,  L and  VWCOS  are  defined  as  before  (I  is  the  from  room, 
J is  the  to  room,  K is  the  vent  number  from  1 to  4,  W is  the  width  of  the  vent,  H is  the 
soffit  or  top  of  the  vent,  L is  the  sill  or  bottom  of  the  vent  and  VWCOS  is  the  wind  velocity 
direction  cosine  ranging  from  1 to  —1).  The  new  parameter  HI  is  the  distance  between  this 
vent  and  the  reference  point  in  room  I.  Similarly,  HJ  is  the  distance  between  this  vent  and 
the  reference  point  in  room  J. 

A .5  Examples 

Figure  12  illustrates  a scenario  involving  a fire  room,  a corridor  and  two  side  rooms.  Room 
1 contains  a constant  100  kW  fire.  The  vents  connecting  the  corridor  (room  2)  to  rooms 
1,  4,  5 and  3 (room  5 is  designated  as  the  outside)  are  located  at  0 m,  1 m,  9 m and  10  m 
respectively.  These  are  specified  as  open  doors  using  the  vent  commands 

HVENT  121. 79  2. 34  0. 00. 0 0.00.0 
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Figure  12:  Diagram  of  a CFAST  scenario  involving  a 10  m corridor  and  three  side  rooms. 

HVENT  2 4 1 .79  2.34  0.0  0.0  1.0  0.0 

HVENT  2 5 1. 79  2. 34  0.0  0.0  9.0  0.0 

HVENT  2 3 1 .79  2.34  0.0  0.0  10.0  0.0 

The  (i,j,k)  vent  coordinate  is  defined  in  Figure  12.  Each  vent  extends  from  the  floor 
(first  0.0)  to  the  ceiling  (2.34)  and  is  0.79  m wide.  The  wind  at  the  vent  is  0.0  m/s  (second 
0.0  on  the  command  line).  Since  room  2 is  a hall,  the  vent  distance  is  the  9’th  parameter 
for  the  first  HVENT  command  and  the  8’th  parameter  for  all  other  HVENT  commands. 

This  scenario  is  simulated  for  30  seconds.  A second  scenario,  the  same  as  the  first 
except  that  the  corridor  is  20  m instead  of  10  m is  also  simulated.  Upper  layer  temperatures 
for  room  3 are  presented  in  Figure  13  along  with  simulations  for  the  same  scenario  modeled 
without  the  hybrid  option.  The  delay  enforced  by  the  hybrid  option  causes  the  delay  in 
temperature  rise  for  room  3 in  the  two  cases. 
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10  m Corridor 


(a)  10  m corridor 


20  m corridor 


(b)  20  m corridor 

Figure  13:  Plot  of  temperature  excess  (°C)  vs.  time  for  corridors  of  length  10  m and  20  m. 
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